Effects of the particle-particle channel on properties of low-lying vibrational states 
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Making use of the finite rank separable approach for the quasiparticle random phase approximation 
enables one to perform nuclear structure calculations in very large two-quasiparticle spaces. The 
approach is extended to take into account the residual particle-particle interaction. The calculations 
are performed by using Skyrme interactions in the particle-hole channel and density-dependent 
zero-range interactions in the particle-particle channel. To illustrate our approach, we study the 
properties of the lowest quadrupole states in the even-even nuclei ^^*Pd, ^^°Cd, ^^''"""Sn, i^s-iserp^ 
and i^^Xe. 
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I. INTRODUCTION 

The low-energy spectrum is a key character of exci- 
tations of nuclei in the presence of pairing correlations. 
The new spectroscopic studies of exotic nuclei stimulate 
a development of the nuclear models [l], 0, Q to de- 
scribe properties of nuclei away from the stability line. 
One of the standard tools for nuclear structure studies is 
the quasiparticle random phase approximation (QRPA) 
with the self-consistent mean-field derived by making use 
of effective nucleon-nucleon interactions which are taken 
whether as non-relativistic two-body forces [1, [H] or de- 
rived from relativistic lagrangians [6[ . Since these QRPA 
calculations are performed with the same energy func- 
tional as that determining the mean-field one does not 
require to introduce new parameters. Such an approach 
describes the properties of the low-lying states less ac- 
curately than more phcnomcnological ones, but the re- 
sults are in a reasonable agreement with experimental 
data B&SKEI. 

When the residual interaction is separable the 
QRPA equations can be easily solved no matter 
how many two-quasiparticle configurations are involved. 
Starting from an effective interaction of the Skyrme 
type, a finite rank separable approximation was pro- 
posed for the particle-hole (p-h) residual interaction. 
Such an approach allows one to perform structure cal- 
culations in very large particle-hole spaces. Thus, the 
self-consistent mean field can be calculated within the 
Hartree-Fock (HF) method with the original Skyrme in- 
teractions whereas the RPA solutions are obtained with 
the finite rank approximation for the p-h matrix ele- 
ments. This approach can be extended to include the 
pairing correlations within the BCS approximation p^ . 
Alternative schemes to factorize the p-h interaction have 
also been considered in [l^, [l^, ■ 

Due to the anharmonicity of vibrations there is a cou- 
pling between one-phonon and more complex states 
and the complexity of calculations beyond the standard 
QRPA increases rapidly with the size of the configura- 
tion space. We have generalized our approach to take 
into account a coupling between the one- and two-phonon 



components of wave functions in Ref. [Tsl , where we fol- 
low the basic ideas of the quasiparticle-phonon model 
(QPM) However, the single-quasiparticle spectrum 
and the parameters of the residual interaction are calcu- 
lated with Skyrme forces. Note that the QPM can 
achieve very detailed predictions for nuclei away from 
closed shells [l^, but it is difficult to extrapolate the 
phcnomcnological parameters of the model hamiltonian 
to new regions of nuclei. 

In the present work, we propose an extension of our 
approach by taking into account the particle-particle (p- 
p) residual interaction. As an application we present 
the evolution of lowest quadrupole states in even-even 
nuclei around the ^^^Sn region. Using the neutron-rich 
radioactive ion beams, the recent B(E2) measurements 
through Coulomb excitation in inverse kinematics give 
an opportunity to compare our results and the experi- 
mental data [23, [m . 

This paper is organized as follows: in Sec. II we sketch 
our method, where the residual interaction is obtained by 
the finite rank approximation. The hamiltonian is con- 
structed in Sec. IIA whereas detailed expressions for the 
residual particle-particle interaction arc in Appendix A. 
We consider the QRPA equations in the case of separa- 
ble residual interactions in Sec. IIB, and the solving of 
these equations is explained in Appendix B. In Sec. Ill 
we show how this approach can be applied to treat the 
low-lying states. Results of calculations for characteris- 
tics of the 2^ states in the Sn, Te isotopes and the N=82 
isotones are discussed in Sec. IV . Conclusions are drawn 
in Sec. V. 



II. METHOD OF CALCULATION 

A. The model hamiltonian 

The starting point of the method is the HF-BCS cal- 
culation 3 of the ground states. We restrict the present 
discussion to the case of spherical symmetry. The contin- 
uous part of the single-particle spectrum is discretized by 
diagonalizing the HF hamiltonian on a harmonic oscilla- 



2 



tor basis [23] ■ We work in the quasiparticlc representa- 
tion defined by the canonical Bogohubov transformation: 

<n = "J<n + (-l)'"''""j«J-m, (1) 

where jm denote the quantum numbers nljm. We use 
the Skyrme interaction [2^ in the p-h channel, while the 
pairing correlations are generated by a surface peaked 
density-dependent zero-range force 

Vpa^r{rl,r2) = Vo (^1 - ^) <5 (n - , (2) 

In definition p{ri) is the particle density in coordi- 
nate space, Pc is equal to the nuclear saturation density, 
the strength Vb is a parameter fixed to reproduce the 
odd-even mass difference of nuclei in the studied region. 
In Sec. Ill, we discuss how to make the choice of the pa- 
rameter Vo- In order to avoid divergences, it is necessary 
to introduce a cut-off in the single-particle space. This 
cut-off limits the active pairing space above the Fermi 
level. As proposed in Ref. [21, [2^, we have used the 
smooth cut-off by multiplying the p-p matrix elements 
with cut-off factors rjj taken as 

ryj =1^1 + exp(^ ) j . (3) 

Ej are the single-particle energies, Xn,p is the chemical 
potential. In our calculations we have set the energy cut- 
off AE equal to 10 MeV above the Fermi level, the width 
parameter /x being 0.5 MeV. 

The residual interaction 1/^^ in the p-h channel and 
yPP, in the p-p channel can be obtained as the second 
derivatives of the energy density functional with respect 
to the particle density p and the pair density^ re- 
spectively. Following the method introduced in [13[ we 
simplify V^^^ by approximating it by its Landau-Migdal 
form. For Skyrme interactions the Landau parameters 
are functions of the coordinate r and all parameters with 



I > 1 vanish. We keep only the I = terms in V^Jl and 

the expressions for Fq^,Gq'^,Fq'^,Gq'^ in terms of the 
Skyrme force parameters can be found in Ref. [2^. In 
this work we study only normal parity states and one 
can neglect the spin-spin terms since they play a minor 
role. The Coulomb and spin-orbit residual interactions 
are also dropped. Therefore we can write the residual 
interaction in the following form: 

C,(ri,r2) = iVo-Mi^o(n) + ^o"(^i)n-T2]<5(ri-r2), (4) 

where a is the channel index a — {ph,pp}, Ti is the isospin 
operator, and TVg = 2kpm* /tt^It? with kp and m* stand- 
ing for the Fermi momentum and nucleon effective mass. 
For the p-p channel the expressions for Fq^ and Fq^^ have 
the following forms: 

FS'ir) = InoVo (l ^) , (5) 

J^rW = n"^M.^' (6) 

As a matter of fact, the definition of the pairing force ([2|) 
involves the energy cut-off of the single-particle space to 
restrict the active pairing space within the mean-field 
approximation. This energy cut-off is still required to 
eliminate the p-p matrix elements of the residual inter- 
action in the case of the subshells that are far from the 
Fermi level. The region of influence of the residual p-p 
interaction is confined to the BCS subspace 

= vlIsA^vimmm, (7) 

where all subshells below the energy cut-off are included. 
The two-body matrix elements 

1^234 = y"(?!)t(ri)02(r2)Kes(ri,r2)03(ri)04(r2)dridr2 

(8) 

can be written as: 



JM 

X (ji ||*''r7|| ja) (j2 ll^'^i^ll m) I""" (J1J2J3J4) , (9) 

^jimij2m2j3m3jirrn = J2^jimij2m2 I JM){j3m3jimi \ JM)ilj,Vj2VhVu 
JM 

^ { n J2 A } + {Jl W^'YxW J3) (J2 IH'nII h) /''P(J1J2J3J4) (10) 



I 

in the p-h and p-p channels, respectively. In the above matrix element of the spherical harmonics Yj^ [27t . 
expressions, J = V2J+ 1, (jiWi-^YjWjz) is the reduced 
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I°'{iii2jzjA) is the radial integral: 



dr 



XUj, (r)uj, (r)iij3 (r)uj^(r)-2-, (11) 

where the radial wave functions u(r) are related to the 
single-particle wave functions: 



0,™(l) = ^^:y,™,(ri,ai). 



ri 



(12) 



We see that the p-li matrix elements are in the sepa- 
rable form in the angular coordinates. The separability 
of the antisymmetrized p-p matrix elements is proved in 
Appendix A. The radial integrals pT|) can be calculated 
accurately by choosing a large enough cut-off radius R 
and using a iV-point integration Gauss formula with ab- 
scissas Tfe and weights Wk- Thus, the residual interaction 
can be reduced to a finite rank separable form: 



1 ^ 



aXfi k—1 rg— ±1 



(13) 



We sum over the proton(p) and neutron(n) indexes and 
the notation {r = {n,p)} is used. A change t ^ — r 
means a change p n. k'J''^'^) (k^pP'''")) are the multipole 
interaction strengths in the p-h (p-p) channel and they 
can be expressed as: 



{ph,k) 



-iRwk 



( Ftivk) \ 
IK'irk) 



(14) 



The multipole operators entering the normal products in 
Eq. (fT3|l are defined as follows: 



M 



(r) = A- 



E 

. ./ / 

J J mm 



(-1) 



j-\-m 



x{jmj - m I Xn)f. 



(Xk) + 

, - a ■ a ' I . 



(15) 



M 



(pp,k) + 
A/1 



(r) = (-l)^-^A-i 



E 



where fj^^] are the single-particle matrix elements of the 
multipole operators: 



{jmj'm \ Xfi)f!j^!"^Tjjrij,a+^at^^,, (16) 



(Xk) 



Uj,{rk)uj^irk){ji\\i^Yx\\j2 



(17) 



The residual interaction (jl3[) is represented in terms of 
bi-fermion quasiparticle operators and their conjugates: 

mm 

(18) 



{jmj m \X^i)aJ^a+,^^,. (19) 



Thus, the hamiltonian of our method has the same form 
as the hamiltonian of the QPM [l2|, but the single- 
quasiparticle spectrum and the parameters of the residual 
interaction are calculated with the Skyrme forces. 



B. QRPA equations for separable residual 
interactions 



We introduce the phonon creation operators 

-(-1)^-^1-^? AC?/; A -A^)), (20) 



where the index A denotes total angular momentum and 
fx is its z-projection in the laboratory system. One as- 
sumes that the ground state is the QRPA phonon vacuum 
I 0). We define the excited states as Q\ ^ \ 0) with the 



normalization condition 

{^\[Qx,^,Ql^.]\0)^5^,'. 



(21) 



Making use of the linearized equation-of-motion ap- 
proach [H one can get the QRPA equations Q: 



A B 

-B -A 



(22) 



where the Ai^"^ -i m ■ / % matrix is related to forward-going 
graphs and the sJ'*'' , s matrix is related to backward- 

^ (jlJl)U2j2)^ 

going graphs. The dimension of the matrices A, B is the 
space size of the two-quasiparticle configurations. In our 
case, we obtain 



( (ph.k 



N 



^ ■'jij'i i2i'2 

k=l 



^ 1 ' Jifi J2J2 



(+) 
32 j 2 



3l3i 3232 



(23) 



B 



(A) 

Ul>j'l)^{j2>j2)qT 



= -x-^ (1 + 5,,,- 



N 



AXk) AXk) , (i 

2^ J 310'J 323'2 ^'^0 



^1 ' 3l3i 3232 



k=l 



{Kq + qi^i )V]i3iV3232 

V 3l3i 3232 3l3i 3232 



(24) 
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where ejj> = Sj+e-j,, r]jy = Vj+V]'^ ^jy = "j^'i' 
and v'"^) = UjUji ± VjVji . The expUcit sohition of the 
corresponding QRPA equations is given in Appendix B. 
Thus, this approach enables one to reduce remarkably 
the dimensions of the matrices that must be inverted to 
perform structure calculations in very large configuration 
spaces. It is shown that the matrix dimensions never 
exceed 6A^ x 6A^ independently of the configuration space 
size. If we omit the residual interaction in the p-p channel 
then the matrix dimension is reduced by a factor 3 (isl . 

III. DETAILS OF CALCULATIONS 



TABLE I: Properties of the 2'^ state in ^^°Te as an illus- 
trative example to demonstrate of effect of the residual p-p 
interaction. 



residual 


E 


B{E2 t) 


interaction 


(MeV) 


(e^ fm*) 


ph 


f.49 


3400 


ph+PP 


l.f5 


4000 


ph-l-pp (cut-off) 


1.27 


3600 



Te, for example. 



We apply our approach to study characteristics of the 
lowest vibrational states in the nuclei around the ^'^^Sn 
region. In this work we use the parametrization SLy4 (23j 
of the Skyrme interaction. One peculiarity is that the pa- 
rameters of the force have been adjusted to describe the 
pure neutron matter. It follows that this parametriza- 
tion is a good candidate to describe isotopic properties 
of nuclei from the /3-stability line to the neutron drip 
line. In our calculations the single-particle continuum is 
discretized (2^ by diagonalizing the HF hamiltonian on 
a basis of twelve harmonic oscillator shells and cutting 
off the single-particle spectra at the energy of 100 MeV. 
This is sufficient to exhaust practically all the energy- 
weighted sum rule within the QRPA. We use the isospin- 
invariant surface-peaked pairing force The value 
0.16fm ~^ is the nuclear saturation density for the SLy4 
force. The pairing strength Vq is fitted to reproduce the 
pairing energies given by 

Pn = \{B {N, Z) + B{N + 2,Z)- 2B {N - 1, Z)) 

(25) 

for neutrons, and similarly for protons. The strength Vq 
is taken equal to -940 MeVfm'' in order to get a rea- 
sonable description of the energies (pS)) for both protons 

and neutrons. The Landau parameters fg'*, Fq^ , Gq^, 

Gn^ expressed in terms of the Skyrme force parameters 
[2q de pen d on kp. As it is pointed out in our previous 
works [Tsl . [T^ one needs to adopt some effective value 
for kp to give an accurate representation of the original 
p-h Skyrme interaction. For the present calculations we 
use the nuclear matter value for kp. Our previous inves- 
tigations [H, [III enable us to conclude that A^=45 for 
the rank of our separable approximation is enough for 
multipolarities A < 6 in nuclei with A < 208. 

It is worth to mention the significance of the energy 
cut-off of the single-particle space to confine the active 
space of the residual p-p interaction. Our choice for the 
cut-off eliminates matrix elements pU)) coupling single- 
particle states inside and outside of the BCS subspace. 
As it is seen from Table I omitting the energy cut-off 
would lead to an overestimation of the effect of the resid- 
ual p-p interaction on the 2^ energy and B{E2 t) in 



IV. RESULTS 
A. Sn isotopes 

As the first application of the method we investigate 
the p-p channel effects on energies and transition prob- 
abilities of 2^ states in ^^^~^^^Sn. Results of our calcu- 
lations for the 2^ energies and B{E2) transition proba- 
bilities are compared with experimental data [13, [2l|, [2^ 
in Fig.l. As it is seen from Fig.l there is a remarkable 
increase of the 2^ energy and B{E2 t) in ^^^Sn in com- 
parison with those in ^^°'^'^^Sn. As it was explained in our 
previous paper [l^ such a behaviour of B{E2 t) is related 
with the proportion between the QRPA amplitudes for 
neutrons and protons in Sn isotopes. Including the p-p 
channel changes contributions of the main configurations 
only slightly, but the general structure of the 2^ remains 
the same. The neutron amplitudes are dominant in all Sn 
isotopes and the contribution of the main neutron config- 
uration {1/111/2,1/111/2} increases from 58% (61% in the 
case of the inclusion the p-h interaction only) in i24Sn to 
85.6% (85.3% for the p-h case) in ^^OSn when neutrons fill 
the subshell lft.ii/2- At the same time the contribution of 
the main proton configuration {2^5/2, 139/2} is decreas- 
ing from 15% in ^^''Sn to 7% in ^^^Sn. The closure of the 
neutron subshell l/iii/2 in ^"^^Sn leads to the vanishing of 
the neutron pairing. The energy of the first neutron two- 
quasiparticle pole {2/7/2, 1/111/2} in ^"^^Sn is larger than 
energies of the first poles in i30,i34gj^ ^^^^^ contribution 
of the {2/7/2, 1/111/2} configuration in the doubly magic 
^■^^Sn is about 61%. Furthermore, the first pole in ^^^Sn 
is closer to the proton poles. This means that the con- 
tribution of the proton two-quasiparticle configurations 
is larger than those in the neighbouring isotopes and as 
a result the main proton configuration {2^5/2,159/2} in 
^■^^Sn exhausts about 33%. In ^'^^Sn the leading contribu- 
tion (about 96%) comes from the neutron configuration 
{2/7/2, 2/7/2} and as consequently the B{E2) value is re- 
duced. Such a behaviour of the 2f energies and B{E2) 
values in the neutron-rich Sn isotopes reflects the shell 
structure in this region. As one can see from Fig.l the 
inclusion of the p-p channel results in a reduction of ener- 
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FIG. 1: Energies and B(E2)-values for up-transitions to the 
first 2+ states in ^^^-"''Sn. 

gies and transition probabilities. The calculations repro- 
duce very well a general behaviour for energies and tran- 
sition probabilities, but there is some overcstimation in 
comparison with experimental data. One can expect an 
improvement if the coupling with the two-phonon com- 
ponents of the wave functions [ll] is taken into account. 
Such calculations are now in progress. It is worth to men- 
tion that the first prediction of the anomalous behaviour 
of 2+ excitations around ^"^^Sn based on the QRPA calcu- 
lations with a separable quadrupole-plus-pairing hamil- 
tonian has been done in [231 ■ Other QRPA calculations 
with Skyrme HO] and Gogny [3l| forces give similar 
results for Sn isotopes. 

B. Te isotopes 

Let us now discuss the Te isotopes. The calculated 2^ 
state energies and transition probabilities B(E2) in the 
i28-i36rpg isotopes and experimental data [10, HU, [H arc 
shown in Fig. 2. The general behaviour of energies of the 
Te isotopes is similar to that of the Sn isotopes. They 
have a maximal value at N=82, but the behaviour of the 
B(E2)-valucs is different and corresponds to a standard 
evolution of the B{E2) near closed shells. As it is seen 
from Fig. 2 there is a decrease of the 2^ energies due to 
the inclusion of the p-p channel. At the same time the 
_B( £'2) -values do not change practically. It means that 
the collectivity of the 2^ states is reduced. The neutron 
configurations exhaust about 17% and 28% of the wave 
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FIG. 2: Energies and B(E2)-values for up-transitions to the 
first 2+ states in i28-i36rpg^ 



function normalization in ^'^^Te and ^^^Te respectively. 
In ^^"^Te the contribution of the neutron configurations 
is less than 3% and the dominant proton configuration 
{157/251.97/2} gives a contribution of about 65% that is 
almost twice larger than in the neighboring Te isotopes. 
As far as a contribution of this configuration into the 
transition probability is less than contributions of other 
proton configurations the i?(£'2)-values have such a be- 
haviour near N=82. The structure of the 2^ in ^'^^Te is 
similar to that in ^^^Te and as a result the i3(£'2)-values 
differ slightly. 

Our calculations describe correctly the isotopic depen- 
dence of energies and transition probabilities and they are 
in a reasonable agreement with the available experimen- 
tal data. It is worth to mention that the anharmonic- 
ity effects are strong for the light Te isotopes and the 
QRPA is not very good in such a case. The i3(£'2)-value 
in the neutron-rich isotope -'^'^^Te is only slightly larger 
than that of ^'^^Te, in contrast to the trend of Ce, Ba 
and Xe isotopes [13, HH, [l^ . Such a behaviour of B{E2) 
is related with the shell structure in this region and an in- 
terplay between the QRPA amplitudes for neutrons and 
protons in Te isotopes. 



C. N=82 isotones 

It is interesting to study a change of the structure of 
the 2'^ states along the N=82 isotones chain. The N=82 
isotones below the doubly magic nucleus ^^^Sn are cru- 
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FIG. 3: Energies and B(E2)-values for up-transitions to the 
first 2^ states in the N=82 isotones. 



TABLE II; {M„/Mp)/{N/Z) ratios for the first 2^ states 



Nucleus 




i3°Cd 


"^Sn "*Te 


i3«Xe 


Theory 


0.47 


0.49 


0.81 0.54 


0.55 



2+ states can be extracted from the proton scattering 
experiments (for example, see [sj]) by looking at the ra- 
tio of the multipolc transition matrix elements Mn/Mp 
that depends on the relative contributions of the pro- 
ton and neutron configurations. Results of our calcula- 
tions are given in Table II, where the M„/Mp ratio for 
i28pd, i3°Cd is less than half of N/Z value, indicating a 
very strong proton contribution. According to our cal- 
culations there is a sharp increase of Mn/Mp at Z=50, 
N=82. Such a behaviour of the multipole transition ma- 
trix elements Mn/Mp in other nuclei can indicate a shell 
closure. Another quantity that characterizes the 2^ state 
is the transition density. As an example the neutron and 
proton transition densities of the 2f state of ^^°Cd are 
displayed in Fig. 4. The neutron transition density is 
shifted outwards as compared to the proton transition 
density due to the presence of the neutron skin. We get 
the similar tendency in the case of the other isotones but 
this effect becomes weak in ^'^^Xe. 




FIG. 4: Neutron (solid line) and proton (dashed line) transi- 
tion densities of the 2^ state of ^^"Cd. 



cial for stellar nucleosynthesis [34|. Results of our calcu- 
lations and existing experimental data [2l], [2^, [H, are 
shown in Fig. 3. It is seen that the inclusion of the p-p 
channel does not change energies and transition prob- 
abilities along this chain. Going along the N=82 iso- 
tones chain one can find that 2'^ states in ^^^Pd and 
^^°Cd have a non collective structure with a domination 
of the proton configuration {lgg/2, ^99/2}- In ^^^Sn as it 
was discussed above the main configurations are the neu- 
tron {2/7/2, 1/111/2} (61%) and the proton {2^5/2, 159/2} 
(33%) ones. In ^^^Te and ^^'^Xe the 2+ states are very 
collective and many proton configurations contribute in 
their structure. The structure pecularities are reflected 
in the B(E2) behaviour in this chain. Higher collectivity 
results in an increase of the transition probability. An 
additional information about the structure of the first 



V. CONCLUSIONS 

A finite rank separable approximation for the QRPA 
calculations with Skyrme interactions that was proposed 
in our previous work is extended to take into account the 
residual particle-particle interaction. This approxima- 
tion enables one to reduce considerably the dimensions 
of the matrices that must be inverted to perform struc- 
ture calculations in very large configuration spaces. As 
an illustration of the method we have studied the ener- 
gies and transition probabilities of the 2f states around 
the ^^^Sn region. Using the same set of parameters we 
describe available experimental data and give predictions 
for the N=82 isotones that are important for stellar nu- 
cleosynthesis. Including the quadrupole p-p interaction 
results in a reduction of the collectivity and this may be 
more important for nuclei far from closed shells. Such cal- 
culations which take into account the two-phonon terms 
in wave functions arc in progress now. 
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APPENDIX A the separable form in the angular coordinates. 

In this appendix, we derive the formulas which help us 
to represent the antisymmetrized p-p matrix elements in In Eq.(lO) the sum over A can be trasformed into: 



2+j3 + J+A ) .74 J 

1 ii h A 



(ji II* ^a|| iz) (h \\i Yx\\ ji) = jij-ihjA (167r) i 



1 Js+U-h-h 



X (l + (-l) 



J3 j4 ^ W jl j2 J 



+ (-1)'^+'*) (-1) 



h+33 / h J k\ i jl J h 



.10- 

2^2 



-1 i 

2^2 



(Al) 



Then, the antisymmetrized p-p matrix elements take the form: 
^12^34 - W243 = ^{jimij2m2 I JA/)(j3m3j4m4 | JM) (ji ||i"'>j|| ^2) (j3 |K^^j|| ^4) {iii2j3ji) VjiVj2VnVj4- 

. m 



JM 



APPENDIX B 

Taking into account the residual p-p interaction we 
show how the finite rank separable form of the residual 
force can simplify the solution of the QRPA equa- 
tions ([22|) . In the 6A^-dimensional space we introduce a 

1^0 (r) 

vector I T>+ (r) I by its components: 
P_ (r) 



n' 



D'f> (r) 

'13 



The index k runs over the TV-dimensional space 
(k=l,2,...,N). Following our previous paper [iJl the 
QRPA equations (22) can be reduced to the set of equa- 
/? = {0,+,-} (Bl) tions: 



/ Moo (r) - 1 Mo+ (r) Mo- (r) 
M+o (r) M++ (r) - 1 M+- (r) 
\ M-o{t) M-+{t) M—{t)-1 




= 0, 



where M is the 2N x 2A'' matrix 



M'^, (r) = 



^0' (r) 



(4 



,l<k,k' <N 



r 



(B2) 



(B3) 



In the definition (|B3]), = = form: 

^(pp,fe )^ matrix elements T'^'^' have the following 



8 



J J 



n' 



One can see that the matrix dimensions never ex- 
ceed 6-/V X 6A^ independently of the size of the two- 
quasiparticlc configuration. The excitation energies ujxi 
are the roots of the secular equation: 



jj' 



To'"' (r) = Tl^o (r) ^ X^f -'^^^ 



where 



Afcfe' _ •'i/ ■'jj' 



/ Xoo (r) - 1 >(o+ (t) TWo- (t) 
det TW+o (r) (r) - 1 X+_ (r) I = 0. 

V M^oir) M-+{t) X__(t)-1 

(B4) 

The phonon amplitudes corresponding to the QRPA 
eigenvalue ujxi arc obtained by Eas.(|22p and the normal- 
ization condition (EU). 
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